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Coaxal Circles and Smith Charts 

The references to Smith charts on page 2 of this report are misleading. The values of z which 
satisfy Equation (11) lie on a circle whose centre is where z = (1 + p 2 )/(l - p) and radius 2p/|l -p 2 |- If 
p - this gives z = 1; as p increases from to 1, the circle, whose centre moves to the right along the real 
axis, increases in radius but remains to the right of the imaginary axis; the circle degenerates into the 
imaginary axis as p -» 1 and for greater values of p it lies to the left of the imaginary axis, and dwindles as p 
increases until as p ~> <x> it reduces to the point z = -1. The length of the tangent from any point z = \k on the 
imaginary axis to any of these circles is (k 2 + l)' 2 whatever the value of p may be, and hence they are correct- 
ly described as coaxal; the radical axis is the imaginary axis. Further, any circle through the points z = ±1 
cuts any of the circles given by Equation (11) for any value of p at right angles. The Smith chart circles are 
a limiting case of the above-mentioned systems when the scale is reduced without limit, so that the circles 
satisfying Equation (11) are reduced to circles touching the imaginary axis at the origin, and the circles 
through the points z = ±1 become circles touching the real axis at the origin. 
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SUMMARY 

One method of broadcast transmitter aerial null-filling consists essen- 
tially in solving an algebraic equation (which in general is of degree 20 or 
more and has complex coefficients) and inverting some of the roots with 
respect to the unit circle in the complex plane. Such inversion does not alter 
the radiation pattern, but may make the practical realization of the aerial 
easier. Hill has also explored possibilities of improving undesirable 
radiation-pattern characteristics by appropriate adjustment of these roots. 
The co-phased quarter-sine amplitude distribution is first discussed in detail; 
this case is specially simple because the above-mentioned algebraic equation 
has real coefficients and a straightforward graphical technique of solution is 
available. Another important practical case is when the roots of the algebraic 
equation are nearly equally spaced round the unit circle in the complex plane. 
The general case is outlined; if the given equation can be regarded as a 
variant of one whose roots are already known (or more easily found) the pro- 
cedure is simplified. 



1. INTRODUCTION 

Various possible methods of aerial null-filling 
are available, and several of them are fully discus- 
sed by Hill 1 . In one of these methods, the equation 

f n (z) = h z"" 1 + l 9 z n - 2 + l 3 z n - 3 + . . . . + /„ = 



(1) 



plays a very important part, where there are n ele- 
ments altogether, and l r is the current in the rth of 
these. In general, the elements will not be co- 
phased, and therefore the coefficients l r must be 
regarded as complex, while n may be in the region 
20 to 30. Thus it appears as if we may be required 
to solve an equation of high degree having complex 
roots. In this note, we first consider the co-phased 
quarter-sine amplitude distribution, which provides 
extremely smooth null-filling but has the disadvan- 
tage that the ratio of maximum to minimum current 
in the tiers is inconveniently large. For this case 



that in which the zero-order approximation to the 
required f n (z) has 



I, = sin i(2r-l)n-/4n} 



(2) 



so that all the I r are real. This means that the 
equation (1) can be solved by a library routine on 
the Research Department Elliott 803B computer for 
an adequate value of n, and in this case a graphical 
method of solution is also available, as fully dis- 
cussed for the particular case n = 8 in Section 2. 
Another case of considerable practical interest is 



k=ln =l,/ 2 = / £ 



/„.! = (3) 



and the required f n (z) is built up by using Hill's 
technique 1 of "perturbing" some of the roots of (1) 
radially; this case is discussed in Section 3. 
Section 4 outlines ways in which the general case 
of solving (1) for complex l r and high n could be 
tackled, with particular attention to the fact that it 
is much easier to adjust the roots of (1) to achieve 
some desirable end and deduce the coefficients than 
vice versa. In particular, it must be noted that for 
the null-filling application, replacing any root of (1) 
by its inverse with respect to the unit circle in'the 
complex z-plane does not alter the radiation pattern 
of the array, but it will of course alter f n (z) and it 
may significantly improve the possibility of realiz- 
ing the radiation pattern to replace several roots of 
f n (z) = by their inverses. If all the roots of 
f n (z) = are replaced by their inverses, we replace 
l r by /„ + i . r and thus do not effectively alter the 
current distribution, but even allowing for this, 
there are 2 n ' x effectively different possibilities of 
changes by inversion to be taken into consideration. 

A routine for solving (1) in the general case 
when n = 5 has been written and fully tested on the 
Research Department computer, but unfortunately 5 
is not a sufficiently high value of n for the null- 
filling application. 



2. THE CO- PHASED QUARTER-SINE AMPLITUDE 
DISTRIBUTION 

All detailed work will be done for the case in 
which n = 8, but the procedure is easily generalized 
for any other value of n. When n - 8, (1) becomes 



z sin — + z sin ^+ z sin- — + 
32 32 32 



13t7 isn- 
't- z sin + sin = 

32 32 

Multiply through (4) by 2j and let 

t _ „jn-/32 



(4) 



(5) 



We have 



z £11 +■— + — - 

\ z z 

z 7 / 1 1 

= -r(l + 7T+ 7TT + 



e 



14 



) 




(6) 



Summing the geometric progressions on both sides 
of (6) 



M2 1 - (^/z 9 ) 1 - (l/^z 8 ) 



l-(£7z) l-(l/<fz) 



(7) 



and when (7) is cleared of fractions and the value 
of £ from (5) is taken into account, we finally 
obtain 

z 9 + z 8 = (z-l) cot — =10-153(z-l) (8) 

32 



Now equation (8) is of degree 9 whereas equation 
(4) is of degree 7; the two extra roots are 

z = £2 = e jir/w. z = £-2 = e -J77/ie (9) 

The remaining roots of (8) are roots of (4) also; 
they were obtained by the Research Department 
computer (in about a minute) as 

0-71674889 ±l-0025036j 

-0-21412418 + 1-3446045J 

J (10) 

-1-1814887 +0-92333381J 

-1-6038395 



so that the loci along which p is constant are the 
coaxal* circles drawn on a Smith chart. Then (8) 
and (11) can be combined to give 

|z| =[10-153p] 1/8 =£ o say (12) 

and if arg z = it can also be shown that 

1-p 2 

cos = k — 2 where k = [£, + (l/£,)]/2 (13) 

1 + p 

and, if <tj = arg(z - 1), ifj = arg(z + 1), 0<cf>, {fj<n 

£o sin £o sin 

tan = ; tan \jj = (14) 

£o cos 8-1 £ cos 8+1 

The condition that z shall be a root of (1) is that 

arg{z 8 (z + l)/(z - 1)1 = 80 - <£ + </r = 2^ltt( m integer) 

(15) 

as well as (12). /z = corresponds to the spurious 
roots (9), while the genuine roots (10), correspond 
to fi = 1, 2, 3, 4 in that order. 

Given p, £ and k ate deducible immediately and 
thus 0, <j) and ^ can be found from (13), (14) and 
(15). As (80 - + xfi) increases steadily with p, 
approximations to the roots (10) can easily be found 
by plotting (80 - <j> + ifr) against p, and improved by 
interpolation. Alternatively, if plotting accuracy is 
sufficient, a value of p corresponding to a circle 
drawn on a Smith chart can be taken, and for this 
value of p the circle (12) can be drawn, with centre 
the origin and radius £ . At the intersection P of 
these two circles (in the upper half plane), 0, and 
ip can all be measured by joining P to the origin 
and to (±1,0), and hence 0, <f> and i(j can be obtained 
and substituted into (15). 

Hence for co-phased quarter-sine amplitude 
distributions the equation (1) takes a particularly 
simple form. This equation can be solved by means 
of a computer library routine for large values of n, 
but it can also be solved in an elementary and 
straightforward manner for any value of n by means 
of the transformation (10) and the graphical equiva- 
lent of using Smith charts. 



3. ROOTS EVENLY DISTRIBUTED ON THE UNIT 
CIRCLE 

Whatever the value of n, there are two cases in 



The following method of attack would permit 
these roots to be found graphically, by elementary 
procedures. Let 



z -1 
z +1 



(ID 



* The locus of points from which equal tangents can be 
drawn to two circles is a straight line called the radical 
axis. A system of circles any two of which have the 
same radical axis is called coaxal. 



which we can solve (1) at sight, namely 

(i)/i = / 2 = / 3 = =/„ 

(ii) h = 1, /„ = R n e ic f> (0<^<2t7; R real and positive) 
While I?, I 3 . . . . I n .i are all zero. 

In case (1) we have 

z"- 1 

/„(z) = ; z = e i2,xn/n il=l2 . . . .(„ - 1) 

z - 1 

so that the roots are uniformly distributed round the 
unit circle, except in so far as the spurious root 
z = 1 is missing. 

In case (ii) the roots are z = Rt^ 2vTT + & /n 
(v = 1,2 . . . . ri) which are uniformly distributed 
around a circle of radius R. 



Unfortunately in both these cases there would 
be deep minima or even actual nulls in the radia- 
tion pattern; this is just what we wish to avoid. In 
case (i), inverting any root with respect to the unit 
circle leaves it unchanged, and the same is nearly 
true in case (ii) if R ~ 1 while if R is appreciably 
different from 1, the ratio of /i to I n will become 
inconveniently large. Hill 1 has found that it is 
useful to alter appropriate roots radially, so that in 
case (ii) for example, the vth root becomes 
kv R expjj(2i>77 + cf>)/n\ where kv is a real constant, 
depending upon v. Hill shows that only one point 
(say P„) of the radiation pattern is significantly 
affected if k v is altered while the other zeros of 
f n (z) are left unchanged; he shows how to calculate 
k v so that an appropriate change in the radiation 
pattern at P v will result. If therefore the k v are 
calculated for points P„ where nulls are obnoxious 
and k v is left unity at places where they are not, 
we have a new set of roots of f n (z) = from which 
the corresponding coefficients of f n (z) can be 
deduced by mere multiplication - a tedious process 
if n is high but relatively simple compared with the 
solution of a high-degree equation with complex 
coefficients. 



4. THE GENERAL CASE IN OUTLINE 

In the general case, our best hope of obtaining 
a solution is to reduce it to the form 



so that the solution of F n (. z ) = is obvious) and 
G n (z) has coefficients which are as small as 
possible. One possibility in the absence of special 
information is to take the coefficients of F„(z) to 
be the real parts of the coefficients of the given 
equation, since the solution of algebraic equations 
with real coefficients by means of a digital com- 
puter is a standard subroutine which is usually 
satisfactory when the coefficients are not delibera- 
tely chosen so that multiple roots will be present. 
The roots of (16) can be obtained explicitly in 
terms of those of F n (z) = by an adaptation of 
Newton's method 2 , and even when |G„(z)| is too 
large for this adaptation to be completely success- 
ful, we can regard the given equation as a variant 
of the equation whose roots are those thus found 
and repeat the process if necessary. 

Iterative processes are likely to work well if 
the given equation can be reduced to the form (16) 
when F n (z) is given by (17) and n is moderate or 
large. For (16) can then be written for iteration in 
the form 



z +a= pe<2ii + l)in/n 
Kfl 



G n (z K ) 



a n 



(n integer) 



(18) 



and we can take z = - a (for all fi) if no better 
approximation is known. Aitken's well-known 8 2 - 
process is widely applicable to processes of this 
kind; if £#, f k + 1 and ^ + 2 are tnree succes- 
sive approximations to (z + a)" arising from the 
iterative process, then a better approximation still 
is 



^K + 1 ~&£k.+ 2 

2 WC + l~ %K ~^K + 2 



(19) 



and the iterative process (18) can then be repeated 
with z + a = (z') 1/fI initially instead of 0. (19) can 
often give a useful result, leading to a true root of 
the given equation, even if the iterative process 
(18), considered in isolation, is divergent, and the 
value of n is only relevant when G n (z K ) has to be 
determined. The biggest difficulty about the S 2 - 
process (19) is that the numerator and denominator 
may both be very small, but as in (19) we are work- 
ing with the nth power of the quantity really re- 
quired, this difficulty is minimized. 



F n (z) = G n (z) 



(16) 



where F n (z) = is an equation whose roots we 
know (either because we have previously had to 
solve a similar but not identical equation, or 
because 



F n (z) = (z + a) n + 



n , an 



(17) 



5. CONCLUSIONS 

Although, in the general case, solving an alge- 
braic equation of high degree with complex coeffi- 



cients is a formidable task, the equations associa- 6. REFERENCES 
ted with aerial null-filling tend to take one of 

several forms for which a relatively simple pro- 1. BBC Research Department Report in preparation, 
cedure is available. In particular, the co-phased 

quarter-sine amplitude distribution case is one 2. The effect of small variations in the coefficients 

which can be handled by elementary graphical pro- of an equation upon its roots. BBC Research 

cedures alone if necessary. Department Report No. F-027, Serial No. 1966/ 

57. 
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